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vectors attached to the vehicle f s marked axes defining the 

spacecraft coordinate system; a*i , ~&2> and ¥3 define the roll, 
pitch, and yaw axes, respectively 

initial orientation of the spacecraft coordinate system 

3 x 3 direction cosine matrix defining the orientation of the 
spacecraft coordinate system IT with respect to the inertial 
coordinate system s’ 

time derivative of matrix A as 

initial value of matrix A as 

3x3 direction cosine matrix defining the orientation of the 
commanded spacecraft coordinate system d with respect to the 
inertial coordinate system s’ 

inner and outer gimbal angles of star tracker 3 

3x3 nonorthogonal matrix defining the orientation of the composite 
tracker coordinate system T with respect to the spacecraft 
coordinate system a 

3x3 nonorthogonal matrix defining the orientation of the composite 
tracker coordinate system T with respect to the inertial 
coordinate system If 

unit eigenvector of the spacecraft-attitude-error matrix 

components of the eigenvector c" in the body roll, pitch, and yaw 
axes , respectively 

cosine 0 

vectors defining the commanded position of the spacecraft coordinate 
system 

vector angular momentum of the system 

moment of inertia of the spacecraft measured about the "ax, and 

113 axes, respectively 

maximum value of I a j (I a3 for this paper) 


X A 11 triads of vectors are orthonormal and form a right-handed set unless 
otherwise stated. 
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I w j moment of inertia of the roll, pitch, and yaw reaction wheels, 

respectively 



gain terms associated with the control law about the body roll, 
pitch, and yaw axes, respectively 


k 2 rate-gain term associated with the control law 


K^,K 0 inner and outer gimbal angles of star tracker 2 

R 3x3 direction cosine matrix defining the orientation of the 

spacecraft coordinate system aT with respect to the commanded 
coordinate system d (error matrix) 

s 0 sine 0 
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large 
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aj 


vectors defining the inertial coordinate system 

composite tracker coordinate system; nonorthogonal triad of vectors 
forming a right-handed system; T2 and T3 are the lines of sight 
to stars 2 and 3, respectively; T\ = T 2 * T3 

vectors defining the orientation of star tracker 2 

vectors defining the orientation of star tracker 3 

control voltage to the motors of the spacecraft reaction wheels 
about the body's marked roll, pitch, and yaw axes 

maximum allowable value of voltage to reaction-wheel motor 

angle of rotation about the desired Euler axis 

magnitude of 0 at which Z- saturates; for this investigation, 

*5 " 0 - 5 ° 

0° <. $ <. 180°J Notation used only with strip chart recorder; the 
> behavior of $ <. 5° is shown with a different 
0° <_ $ <. 5° ) recorder gain than 0° <. 0 <_ 180° 

yaw, roll, and pitch Euler angles defining the orientation of the 
air-bearing table; order of rotation is as shown 

inertial angular rates of the body about the ai , a*2, and i T 3 axes, 
respectively 


w am maximum component of body inertial rate 


w c j computed values of oj a j ; it refers to rates computed from the 

attitude matrix or from tachometers 
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roll, pitch, and yaw reaction-wheel angular rates, respectively 
initial speeds of three reaction wheels 

Subscript 

axis system components, j = 1, 2, 3 



ANALYSIS AND VERIFICATION OF A SPACECRAFT 


WIDE-ANGLE ATTITUDE -CONTROL SYSTEM 
Francis J. Moran and Bruce H. Dishman 
Ames Research Center 


SUMMARY 


Attitude control to ensure stability of a spacecraft sluing through 
large angles was experimentally investigated. The investigation examined a 
control law synthesized from Euler 1 s theorem on rotation, which states that 
any attitude change of a rigid body may be accomplished by a simple rotation 
about an appropriate axis. The theoretical analysis showed that a system 
using the control law is globally stable. 

The purpose of this investigation was to verify by means of an air- 
bearing table simulation the behavior of a system using the control law. It 
also was desired to compare the operation of the control law when vehicle rate 
information was obtained from (1) rate gyroscopes, (2) the spacecraft-attitude 
direction cosine matrix, and (3) tachometers mounted on the reaction wheels. 
The system using rate-gyro information was tested with and without initial 
angular momentum; the other two methods were tested only with zero initial 
momentum. 

The experimental apparatus used in this analysis was an air-bearing 
table representing the spacecraft, three reaction wheels for control, three 
rate gyros for vehicle rate information, and two gimbaled star trackers for 
position information. An online digital computer performed as a real-time 
controller. The computer calculated the desired axis of rotation from the 
star-tracker information and generated simultaneous commands to the three 
torquers. The three control signals were converted to analog form and 
transmitted to the reaction-wheel torquers. 

The results of the simulation show that the system using the control law 
performed as predicted by theory. The air-bearing table rotated about the 
desired axis while sluing through angles of 55° in yaw and 10° in pitch and 
roll. The combined steady-state error about all three axes was 0.1°, which 
represents the sensitivity limit of the attitude sensors. 


INTRODUCTION 


The need for a system that ensures stability while controlling the 
rotational position of a rigid body sluing through large angles is encoun- 
tered in a program such as the Orbiting Astronomical Observatory (0A0) . The 
problem is to point a telescope, rigidly attached to the spacecraft, in any 
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direction dictated by a ground-based observer. The present 0A0 program 
accomplishes wide-angle sluing by commanding the vehicle to complete a 
sequence of open-loop, single-axis maneuvers (ref. 1). The slue-control sys- 
tem reorients the vehicle to within the limits of the coarse control system. 
The coarse control system then stabilizes the satellite by resolving the star- 
tracker gimbal angles into the proper components to drive the vehicle so that 
the target star is within the field of view of the telescopes. Reference 2 
describes another attitude-control system that uses the star-tracker gimbal- 
angle information to stabilize the satellite. 

The purpose of this investigation was to simulate and verify a control 
law that could be used for both stabilization and reorientation. The control 
law was synthesized from Euler’s theorem on rotation, which states that any 
attitude change of a rigid body may be accomplished by a single rotation about 
an appropriate axis. It also was desired to compare the behavior of the sys r 
tern when the rate information required by the control law was obtained from 
gyros, the attitude direction cosine matrix, and tachometers mounted on the 
reaction wheels. 

The control-system simulation was performed on a three-degrees-of- 
freedom air-bearing table. The experimental apparatus used with the table 
consisted of two gimbaled star trackers as position sensors, reaction wheels 
as torquers, rate gyros for damping, and a digital computer as an online 
controller. 

The report first describes the experimental apparatus used to verify the 
control law. The attitude-error matrix equation is then developed followed by 
a brief review of the attitude control law. Then the methods of deriving 
inertial angular body rates from the vehicle-attitude matrix and the reaction- 
wheel tachometers are developed. Finally, the test procedures and results 
describing the performance of a system using the control law are presented. 


EXPERIMENTAL APPARATUS 


The attitude-control system was mechanized on the Ames satellite 
attitude-control simulator (fig. 1), which includes an air-bearing table in a 
vacuum chamber that can be evacuated to a pressure of 1 mm Hg. (The chamber 
was never evacuated during the simulation.) Both the chamber and the air- 
bearing table are mounted on a seismic foundation to isolate the table from 
ground disturbances. The air-bearing table configuration and the control- 
system package used for this analysis is shown in figure 2. The complete 
package weighed 1800 lb; its inertia was 180 slug-ft^ about the yaw axis and 
145 slug-ft 2 about the pitch and roll axes. The air-bearing table has 
unlimited freedom of motion about the yaw axis and is limited to ±29° about 
the pitch and roll axes. (The simulator is described in greater detail in the 
appendix. ) 

The attitude of the table is controlled by three inertia wheels located 
along the control axes of the table (fig. 2). The inertia wheels are driven 
by dc motors with a stall torque of 15 in.-oz. The motor and reaction wheel 
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combination has a time constant of 35 sec and a maximum speed of 3000 rpm. 

The inertia of the wheels was 0.014 slug-ft 2 . The angular velocity of the 
reaction wheel was available from a tachometer. 

The attitude information used to generate the control signals was derived 
from two gimbaled star trackers, which are positioned on the table so that 
their outer gimbal axes are parallel and aligned with the table f s roll axis. 
The trackers have angular freedom of ±50° about both the inner and outer 
gimbal axes. The angular position of the gimbals was determined by measuring 
the output of a resolver, which has an accuracy of ±5.0 min of arc. The star 
trackers have a maximum sluing rate of 0.75° per sec. 

The trackers point at simulated stars located within the air-bearing- 
table chamber. The star simulators have a clear aperture of 12 in. and pro- 
vide a star beam collimated to within ±5.0 arc sec relative to the optical 
axis. The star can be adjusted from -2.0 to +6.5 apparent star magnitude. 

The output from the two star trackers, the three rate gyros, and the 
three tachometers on the table was sent in analog form to a digital computer. 
This computer acts as a real-time controller for the air-bearing table in that 
it accepts the basic analog signals, converts them to digital data, and per- 
forms the necessary matrix operations to generate the specified control law. 
The use of Fortran, facilitated computer changes in the control law. The sys- 
tem program was cycled through the computer at a frequency of about 30 Hz. 

The control signals generated by the computer were transmitted in analog form 
to the reaction wheels. 


DEVELOPMENT OF THE ATTITUDE- ERROR MATRIX EQUATION 


The position information required by the control law makes it necessary 
to derive a measure of the error in attitude. The first section describes the 
development of the attitude error matrix. The second section describes the 
determination of the satellite’s attitude from the star tracker information 
which is necessary to calculate the error matrix. 


Development of the Error Matrix 

The commanded orientation of the spacecraft with respect to the inertial 
coordinate system (fig. 3) is given by:^ 



2 Figure 3 shows the star-tracker inner-gimbal angles at zero. Thus, 
vectors_ T 2 i and T$i are oriented along the two outer gimbal axes. Vectors 
d -2 and d^ are also shown lying in the and a" 2 plane. This arrangement 

is for convenience in drawing and does not represent the general case. 
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The present position of the vehicle with respect to the inertial 
coordinate system is given by: 



From equations (1) 
matrix, it follows that 
the desired position is 


Equation (3) shows that 
matrix product A as A j s 

the error matrix R: 


When the vehicle "a is 
the deviation of R from I is a measure of the error in the attitude of the 
spacecraft. The control law information was extracted from this matrix. 


Determination of Spacecraft Attitude 

The attitude matrix A as must be determined as a basis for calculating 
R; the matrix A ds is given. The attitude of the spacecraft can be defined 
with respect to the lines of sight to two stars. If it is assumed that the 
tracker perfectly tracks the star, then the lines of sight to the stars coin- 
cide with the tracker's optical axes. The attitude of the spacecraft can now 
be defined with respect to the optical axes of the two trackers. For example, 
the relationships between the star-tracker coordinate system and the 
spacecraft are 


and (2) and the fact that A ds is an orthogonal 
the actual position of the vehicle with respect to 
given by: 


= A 




as ds 


(3) 


the vehicle attains the desired attitude when the 
forms an identity matrix. This product is defined as 


R = A as A ds (4) 

in its desired orientation d, then R = I. Therefore, 


and 



X(Ki)V(K 0 )Y(l) 




W(Ai)V(A 0 )Y(2) 



(5) 


( 6 ) 
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where 


0 

0 

1 


X(Ki) 



sK t 

£% 

0 


V(K 0 ) = 



Y(l) = 




W(Ai) 


V (A 0 ) 


Y(2) 




Let T 2 and T3 be unit vectors in the direction of the lines of sight 
of stars 2 and 3, respectively. Then 


t 2 = T 22 
T 3 = T 33 ; 


(7) 


A third vector is defined as 


T x = T 2 x T 3 (8) 

where Tj is not, in general, a unit vector. The Tj., T 2 , T3 coordinate 
system, called the tracker coordinate system, usually is nonorthogonal . The 
relationship between the tracker coordinate system and the spacecraft 
coordinate system is 



In a real application, the location of the stars with respect to inertial 
space would be obtained from some astronomical data source. The location of 
the simulated stars with respect to the inertial space used for the simulator 
is determined from the star-tracker gimbal angles and the table Euler angle 
measurements. For example. 


and 


where 



X(K i )V(K 0 )YCl)A as (K) 




W(Ai)V(A 0 )Y(2)A as (K) 



AasOO 


( c 0 K 0 

0 1 
S0K 0 




s^K 

c^K 

0 



( 10 ) 


( 11 ) 


( 12 ) 


Equation (12) is calculated for known table Euler angles and is used in equa- 
tions (10) and (11) in conjunction with known tracker gimbal angles. Since 
T 2 = T 22 and T 3 = T 3 3 , the position of the tracker coordinate system relative 
to inertial space is 



(13) 


The tracker gimbal angle and table Euler angle measurements are then 
used to calculate B s . The B s matrix locates the simulated stars with 
respect to the inertial frame. If the right side of equations (9) and (13) 
are equated, one obtains 


Comparing equations (14) 


Then the error matrix R, 



and ( 2 ) , one can write 

A as = (B s _ 1 B a ) T 
defined in equation (4) , 
R = (B s _ 1 B a ) T A^ s 


(14) 


(15) 

can be written as 

(16) 
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In a space application, the matrix B s , which defines the orientation of 
the composite tracker coordinate system when locked on the selected guide 
stars to the inertial coordinate system, would be obtained directly from 
as t r onomi cal da t a . 


THE ATTITUDE CONTROL LAW 


The control law was synthesized from Euler’s theorem on rotation (ref. 3), 
which states that any attitude change of a rigid body can be accomplished by a 
single rotation about an appropriate axis. The axis of rotation is defined by 
the eigenvector c* of the error matrix R. The angle of rotation about ' c 
to make the body axes coincident with the desired axis is defined as $. The 
command to the reaction-wheel motors to rotate the spacecraft through an angle 
$ about the eigenvector c* is 

Zj = -kjjfsatO, $ s ) c j + k 2 w a;) .] j = 1, 2, 3 (17) 

where Zj is the control voltage to the jth reaction wheel, co a j is the 
body angular rate about the jth marked body axis, and k^j and k 2 are con- 
stants. The first term in the brackets represents the position information 
and is obtained from the error matrix R determined in the previous section. 
The term sat(0, $ s ) is 


sat($, $ s ) 



for $ >_ $ s 
for $ < $ s 


where the magnitude of rotation $ is 

$ = COS -1 y (Rn + R 22 + R 33 - 1) 


(18) 


and 

*s = °- 5 ° 

The direction of rotation is given by the components of the eigenvector c - 
which are 




C2 


l_ R 23 ~ R 32 
2 sin $ 

L R 31 ~ R 13 
2 sin $ 
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1 R 12 " R 21 
3 2 sin $ 


The Rjj terms refer to the off-diagonal components of the error matrix R. 

The control law can be approximated when the error in attitude is small. 
The approximation is desired because the eigenvector orientation is not 
defined when 4> = 0. Thus, the components of the eigenvector become noisy 


then be written as 


where q- are 


-k 


ij 


5i 


= 


sin 

The 

k 2“aj) 

j 

( R 23 " 

r 32) 

( R 31 - 

R 1 3 ) 

( R 1 2 - 

R 2 1 ) 


j = 1, 2, 3 


(19) 


The control law given by equation (19) is independent of the components of the 
eigenvector. A block diagram of the control law computation is shown in 
figure 4. 

The rate information required by the control law (eqs. (17) and (19)) was 
obtained by three different methods: (1) directly from the rate gyroscopes, 

considered the basic mode of operation; (2) from A a s matrix computations; 
and (3) from the reaction-wheel tachometer measurements. Methods (2) and (3) 
are considered backup systems to the gyroscopes. 

The body rate is derived from the A as matrix according to the 
relationship 


^as ~ s ( w a)A a s 


( 20 ) 


where 


s(o) a ) = 


and A as may be evaluated from 



w a3 

0 

-“al 



( 21 ) 


Aas( n + 1) - Aas (n) 


'as 


( 22 ) 
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which represents the discrete form of the time derivative of Aas; A as (n) 
and Aas (n + 1) ate values of the nine components of the Aas matrix at the 
beginning and end of a time interval specified by T. The matrix equation 
then takes the form: 


If 



T 

B = A as (n + l)A as (n)/T, then the 


w ci = 
U C2 = 
W C3 = 


as (n + 1) - A as (n)“| T 
T A as 

components of the body rates are 

B23 - B32 ' 

2 

B 31 - B 13 
2 > 

Bl 2 - B 21 

2 / 


(23) 


(24) 


where the notation to c , substituted for u> a , indicates a computed value. The 
rates are obtained by summing and averaging the two solutions available for 
each rate component from equation (23) . 

To compute the body rates from reaction-wheel-tachometer information, two 
assumptions must be made: (1) the system is initially stabilized so that all 

the system momentum is stored in the reaction wheels; and (2) the system is 
conservative - that is, the external torques are zero. The total system 
momentum at any time is 

H = ( I ai u) ai + I \n a wO'*i + C I a 2 w a 2 + I W 2 a) w 2^2 + (I a3 “ a 3 + I w3 u w3^3 ( 25 ) 
Since the total initial system momentum is stored in the motors, one can write 

H = Iwi^lCODaxCO) + ^ 2^ 2 (0) + ^W3***W3 (®) ^3 CO) (26) 

The starting position of the body a^(0) can be expressed in terms of the body 
position at any time a^ as 

/ai(0)\ /sA /aA 

( a 2 (0) J = A as (0) f s 2 1 = A as (0) A as (t) ( a 2 j (27) 

\a 3 (0)/ \s 3 / \a 3/ 

Substituting equation (27) into equation (26) and then equating equations (25) 
and (26) , one obtains 
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MP 




where M and N are 1x3 row ratrices whose components are given as 


" ^wj^wj (0) 

N j = I aj aJ aj + Vj^wj 


j = 1 , 2 , 3 


and 


P = A as (0)A^ ft) 


Solving for 
yields 


03 a jj > and denoting it by w c j to indicate a computed value, 
3 

- I wj a) wj 

0) c j = — j-j j = 1, 2, 3 


(28) 


TEST PROCEDURES AND RESULTS 


The theoretical study (ref. 3) shows that a vehicle sluing under the 
influence of the previously given control law will move about the eigenvector 
of the R matrix. The time history of the angle <f> and the components of 
the eigenvector cT are a measure of the control-system performance. If the 
total system momentum remains zero, the time history of the ~c components 
should be a constant value and the angle $ should be a decreasing function 
of constant slope. The final value of $ gives the system-pointing error. 

The results of the simulation are summarized in figures 5 through 11. 

The results were output from the digital computer, through the digital analog 
converters to a strip chart recorder. Data were obtained once every computa- 
tion cycle (35 msec) . A line printer coupled directly to the digital computer 
provided numerical values for the data. A printout was made at the beginning 
and end of each slue. The numerical values are shown in parenthesis on the 
figures. The parameters displayed in the figures include the angle $, the 
eigenvector components c j , the control voltages Zj, and the body rates o) a j • 
Since system-pointing accuracy as well as the wide-angle sluing character- 
istics were of interest, the level of control saturation was set at 
$ s = 0.5°. The magnitude of k X j was determined for each axis from the 
relationship 


k 


l 3 ~ 


^m^aj 

*am 


j = 1, 2, 3 


(29) 
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The gain factor k 2 , which is the same for all three axes, is 



where was set equal to 0.4° per sec; that is, k 2 = 2.5. 


Rate Gyroscopes 

Figure 5 shows the response of the system when rate gyroscopes are used 
for damping. The table is positioned initially at Euler angles of -4.78° in 
roll, -5.74° in pitch, and 74.28° in yaw, and is commanded to slue to +5.0°, 
+5.0°, +20.0°. The direction between the eigenvector and the respective body 
axes is given by the inverse cosine of each c" component. For this run, the 
initial orientation of the eigenvector with respect to the roll, pitch, and 
yaw axes is 81°, 79°55 l , and 13°25 ! , respectively. As expected, the major 
component of the command is about the yaw axis. The time history of the 
eigenvector components shows that they remain nearly constant during the slue 
until $ = 0.5°. When $ j< $ s > t ^ ie approximate control law, which is indepen- 
dent of ~£, is used, and the eigenvector components vary rapidly as indicated 
by their traces in the figure. 

The initial and commanded angles result in an angle $ = 55.94°. The 
time history of the angle $ indicates that it decreases with a constant 
negative slope as predicted by the theory. The recorder sensitivity of $ is 
increased when $ <. 5.0° to show the performance of the system near the target 
point. (The computer was programmed such that $ small was zero until 
$ = 5.0°.) The trace of $ shows that the system has two overshoots and 
settles to within 0.12° of the commanded value. As indicated by the Euler 
angles, the corresponding pointing error about each control axis is (0.02°, 
0.08°, 0.02°). 

From the time histories of the control voltages Zj and the body rates 
co a j , it is observed that the voltage about the axis with a maximum error 

initially saturates and then decreases to a level necessary to maintain the 
maximum desired rate (tog^) . In this case, the yaw axis has the maximum 

error. The average control voltages about the roll and pitch axes are 
proportionally lower as determined by the values of the c“ components. 

The behavior of the system using gyro damping to a sequence of commands 
is shown in figure 6. This sequence is mechanized by commanding the system to 
one position, letting it stabilize at this position, then commanding it to the 
second position, etc. The sequence of commands shows the behavior of the sys- 
tem in a scanning- type operation. Initially, the table's angular position is 
(0.34°, 0.39°, 54.60°) and is commanded to (-5°, -5°, +20°). For this command, 
the major control torque is about the yaw axis. The behavior of Cj and $ is 

quite similar to the run discussed in figure 5. There is no momentum buildup 
during this phase of the run. 


The second sequence commands the table to (+5°, +5°, +20°). During this 
maneuver, the control torques are entirely about the pitch and yaw axes, and 


11 



the eigenvector components vary slightly as a result of undesirable torques 
causing increased system momentum. At t = 178 sec the system is commanded to 
(-10°, -5°, +54°). During this maneuver the c" components vary considerably 
because of the buildup of momentum. The initial system momentum and undesir- 
able external torques cause the eigenvector orientation to change 17° between 
t = 183 sec and t = 260 sec. Nevertheless, the angle approaches zero in a 
nearly linear manner, and the system stabilizes with an error of 0.10°. The 
final command is (0°, 0°, 54°); the system stabilizes to within 0.10° of this 
position as indicated by the final value of $. 


Rate From the A as Matrix 

Figures 7 and 8 show system response when rate damping information is 
derived from the Aas matrix. The procedure used is essentially a time 
differentiation that tends to amplify the system noise to an extent, in fact, 
that it completely masks the rate signal. To reduce the noise, the A as 
matrix is averaged over 100 msec. The averaging process acts as a filter but 
does not introduce enough lag to affect the stability of the system. A com- 
parison of the derived rate a> c j and the rate from the gyroscope to a j shows 
that the derived rate, even when averaged, is more noisy. 

The eigenvector components of the system using the derived rate show more 
variation than those of the system using the gyros (figs. 7 and 8). For 
example, the angle between the eigenvector at t = 0 and t = 125 sec was com- 
puted as 20°. However, both the transient and steady-state responses of the 
system are still quite acceptable. 


Rate From the Reaction-Wheel Tachometers 

Figure 9 shows control system damping in which the body rate information 
is computed from dc tachometers attached to the reaction wheels. Computations 
assumed zero external torques. However, the system is subjected to distur- 
bance torques that cause variations in system momentum. At the initial and 
final positions, the system momentum is in the reaction wheels: the initial 

wheel speeds are zero, but the final wheel speeds are nonzero. The final body 
rate computed from the tachometer measurements, therefore, is nonzero and 
disagrees with the gyro information. The error in the rate signal can be 
observed by comparing the body rate signals derived from the tachometers to 
the gyro measurements on figure 9. At time t = 120 sec, the pitch-rate 
signal computed from the tachometer is of opposite sign from the true body- 
rate measurements. The final rate error is mainly in the pitch axis. From 
the data obtained, this system clearly has the worst response of the three 
tested, because the initial momentum vector changes as the system is subjected 
to disturbance torques. Since the computation is based on the system's ini- 
tial momentum, any momentum change causes an error in the computed rate signal. 
The rate signal error caused the system to develop an additional position 
error as indicated by the $ trace. 

The increased steady-state error can be eliminated through updating the 
system momentum at the end of the maneuver. For example, assume the system 
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begins the maneuver with zero momentum and stabilizes with nonzero momentum. 
After the system stabilizes, one can calculate, by means of equation (26), 
the system momentum. The derived body rate can then be calculated in the 
prescribed manner. The effect of updating the system momentum is to eliminate 
the error in the derived body rate caused by the external torques. The 
results of this procedure are shown in figure 10. The new system momentum and 
the derived rates are recalculated at t = 217 sec. Note that the error in 
the derived body rates is reduced, and that the pointing error is reduced 
from 0.38° to 0.15°. 


Nonzero Initial Momentum With Rate Gyroscopes 

Figure 11 shows the effect of nonzero initial momentum on the behavior of 
the control system using rate gyros for damping. The tachometer traces show 
that the initial rates of the roll, pitch, and yaw reaction wheels were 1600, 
1150, and 1050 rpm, respectively. The final wheel velocities were very close 
to these initial values, which implies a low-magnitude disturbance torque. 
During the run the c. vary as would be expected since the momentum vector of 
the system is varying? The traces of the motor control voltages show the off- 
set necessary to maintain the wheels spinning at the initial rate. Therefore, 
the system performance does not deteriorate significantly when the initial 
system momentum is nonzero. Furthermore the system is well behaved and has a 
pointing error of 0.12°. 


CONCLUDING REMARKS 


The results of the simulation show that the wide-angle attitude- control 
system based on Euler’s theorem on rotation performs as predicted by the 
theoretical analysis. The angular motion is about the Euler rotation axis 
rather than about three axes separately. The system has proven global stabil- 
ity, and the angular rotation of the vehicle is well behaved and predictable 
during a slue. 

The best method for obtaining rate information for the control law is the 
rate gyroscope. External torques or initial momentum that causes a variation 
in the eigenvector during the run does not cause a significant decrease in 
system performance. Inertial body rates also can be determined indirectly, 
which obviates a flight abort in case of rate gyroscope failure. The inertial 
body rates could be computed from either the reaction-wheel tachometer measure- 
ments or from the spacecraft- attitude matrix (A as ) . The rates computed from 
the tachometer measurements are less noisy than the rates obtained from the 
attitude matrix. However, for the tachometer system, a variation in system 
momentum causes an increase in the steady-state pointing error, which can be 
decreased by updating the system momentum vector when the system reaches 
steady state. The performance of the system using the rates derived from the 
attitude matrix is not affected by varying system momentum. 
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The type of control system shown in this mechanization could be used 
very effectively on a vehicle such as the 0A0, where the ability to command 
the spacecraft through large angles is desired. 


Ames Research Center 

National Aeronautics and Space Administration 

Moffett Field, California, 94035, November 25, 1969 
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APPENDIX 


DESCRIPTION OF SATELLITE ATTITUDE- CONTROL SIMULATOR 


The attitude-control system of this report was mechanized on the Ames 
Satellite Attitude-Control Simulator (fig. 1). The air-bearing table is con- 
tained in a spherical chamber, which can be evacuated to a pressure of 1 mm Hg. 
The chamber and air-bearing-table assembly are supported on a seismic 
foundation. 

The table rests on a 7-in. -diameter stainless steel ball supported by a 
cushion of air injected into a mating socket. The gap between the ball and 
the socket is 150 yin. Air is injected into the socket at a pressure of 
70 psia, and is bled off from the socket and routed to control jets located on 
the three axes of the table. With these jets and the control console associ- 
ated with the table, it is possible to control the position of the table inde- 
pendent of the control system package under test. It is possible to add 
inertia cones to the basic table and raise the inertia about the yaw axis to 
750 slug- ft 2 and about the pitch and roll axes to 400 slug-ft 2 . 

Communications with the table and power to the table are transmitted via 
a mercury pool and slip ring assembly. From the table, signals pass through 
individual pins arranged in a circle and located on the table. The pins pro- 
ject into individual pools of mercury located on a following gimbal system. 

(The following gimbal system tracks the table by utilizing a three-axis, auto- 
collimating error-detection system.) From the base of the mercury pools, the 
signals then pass through a slip-ring assembly, located on the yaw gimbal, 
which allows unlimited yaw motion. This technique of getting data on and off 
the table eliminates the disturbance torques that would be caused by hanging 
wires from above the table and the complexity of telemetering the data. 

Also connected to the following gimbal system are position and rate 
encoders that give the table Euler angles and rates in a natural binary form 
at the control console. This information is then converted into a visual 
readout . 
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Figure 2.- Air bearing table with associated equipment. 
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Figure 3.- Coordinate systems orientation. 
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Figure 4.- Block diagram of the control system. 
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Figure 6.- Gyro rate damping - multiple command. 
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Figure 7.- Rate damping derived from A as - single command - Concluded. 
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Figure 9.- Tachometer rate damping - single command. 
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Figure 10.- Tachometer rate damping - momentum buildup compensation. 
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